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In this proceedings we summarize our calculation of the phase diagram of QCD at non-zero 
temperature and chemical potential on S l x S 3 from one-loop perturbation theory [1], which is 
valid in the limit R « Aq^ d , where R is the radius of S 3 . We calculate several observables 
including the Polyakov lines and the quark number, for large number of colors N and large number 
of quark flavors Nf, onS 1 x5 3 , and compare with results for the same system with N = 3, and with 
results for N = 2 lattice QCD. For > 2 the action is complex and the dominant contributions to 
the path integral occur in the space of complexified gauge field configurations. This results in the 
expectation values of the eigenvalues of the Polyakov line lying off the unit circle and out in the 
complex plane. This is an important issue for the lattice, and also for the calculation on S l x 5 3 
in the large limit where we obtain analytical results using the saddle point approximation. It is 
thus necessary to adapt available techniques to locate the stationary solutions in the complexified 
gauge field configuration space. 
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1. Introduction 

Calculation of the phase diagram of QCD as a function of temperature and chemical potential 
is complicated by what is known as the sign problem: the action of QCD becomes complex in the 
presence of a non-zero chemical potential. This prevents conventional methods of lattice simula- 
tion because the Boltzmann factor can no longer take a probability interpretation. If the chemical 
potential is not too large several useful alternative methods of calculating the phase diagram from 
lattice simulations have been developed (for a recent review see [2]). These are valid for \ijT < 1. 
Calculation of the phase diagram at asymptotically large chemical potentials and low temperatures 
is also possible from perturbation theory [3]. The region of the phase diagram that is most difficult 
to access is that for moderate chemical potentials and low temperatures, where the sign problem is 
thought to be more severe, and both conventional lattice techniques and conventional perturbation 
theory are not available (However, simulations using complex Langevin techniques have the poten- 
tial to probe this region. See [4] for a recent report on the progress and issues using this technique 
in the XY model.). 

At the moment, to calculate in the regime of moderate chemical potentials and low tempera- 
tures a sacrifice is necessary. We have opted to sacrifice the large volume limit, compactifying the 
spatial volume such that perturbation theory becomes valid (R << Ag^ D where R is the radius of 
the S 3 ). On a hyperspherical manifold, S l x S 3 (as opposed to (S 1 ) 4 ), small volume results for the 
phase diagrams of related theories, such as Yang-Mills theory and adjoint QCD, qualitatively re- 
semble lattice results (see for example [5,6]). A thermodynamic limit of sorts is obtained by taking 
N — > °o. In this case sharp, well-defined phase transitions are possible, even in finite spatial vol- 
ume. Thus it sometimes happens that larger volume lattice results of a small N theory more closely 
resemble the small volume large N theory, than the small volume theory of the same N. It is useful 
to consider both perturbative results and lattice results to distinguish between small volume effects, 
lattice artifacts, large N effects, and non-perturbative contributions. This proceedings reviews our 
perturbative results for QCD on S l x S 3 [1], comparing with lattice results for N = 2 [7]. 



2. Background 

For our perturbative calculations, all quantities are derived from the 1-loop action of QCD 
on S l x S 3 , which was originally derived in [5, 8] for theories with more general matter content, 
and summarized for QCD in [1]. For QCD with N colors and Nf quark flavors with mass m and 
chemical potential ji and at a temperature T = 1 /j3 the action is given by 

°° i N 

S(e i ) = J £-(l-z v (nP/R)) £ cos(»(ft-0 y )) 

(2.1) 

+ £ ±-tLN f Zf{np/R,mR) £ e npn+ine t + e -nPn-ine i 

where the 0, are the angles of the Polyakov line matrix P = diagje' 9 ' , ...,e' dN }, and z v , Zf are the 
single particle partition functions for (vector) gauge fields and fermions, respectively, defined by 
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Zv (P/R) = Z d f ] ^ et =2^e(l + 2)e-^y R , (2.2) 
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Figure 1: Quark number (Left) and Polyakov lines (Right) as a function of the chemical potential for QCD 
on S l x S 3 . (Right). N = 3, N f = 1, m = 0, P/R = 30 (low T). 



z f (P/R,mR) = £ dfe-^i = 2 £ + i) e -M'+*)W*7* . 



(2.3) 



3. N = 3 results: Quark number and Polyakov lines ^ and ^ i 



We are interested in calculating quantities that are derivable from the partition function, which 
in the low temperature limit takes the form 



Z(P/R) = J [dfl]exp 
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£ - Tr A (P n ) + (-\yN f Zf(nl5/R,mR)e n ^Tr F (P n ) 



(3.1) 



where [dd] = Y\!i = \ddj. We consider a few observables in this proceedings (more are calculated 
in [1,7]), specifically the quark number and the Polyakov lines given by 
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(3.2) 
(3.3) 



For Af = 3 we are able to calculate these observables numerically because there are only the integrals 
over d\ and B 2 to compute since P G SU (3). Figure |l| shows the quark number and Polyakov lines 
as a function of ^R in the massless limit. The level structure in Figure [j] (Left) of the quark number 



is caused by the Fermi-Dirac distribution function in eq. ( |3.2| ) which results from taking the /I 
derivative. In the low temperature limit (j8 — > °°) the Fermi-Dirac distribution function is zero for 
/I < = (/ + A) /R, and 1 for /I > , so as /I passes each energy level the quark number rises 
another step. Thus, the quark number of the Lth level is given by j¥l = AW/Xf =1 2£(£ + 1) . 

Close inspection of Figure [j] (Right) shows that the Polyakov lines have the property 2?\ ^ 
x . This results at non-zero jj. for N > 2 from the fact that the dominant contributions to the 
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path integral lie in the space where the Polyakov line angles 0, are complex, a consequence of the 
complex action. This will be important to our large N analysis since there the saddle point method 
is valid. It is also interesting that the Polyakov lines as a function of p, have a deconfinement spike 

at each energy level Ef \ which results due to the finite separation of the E^f' on a finite volume 

if) 

and indicates that deconfinement only occurs during quark creation. As p is increased past Ef the 
newly created quarks combine and form color singlets. 



4. Large N theory 



Performing the sums over n in eq. (2. 1) we obtain for the action 
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+ £ 21(1 + 2) log [cosh(j8e / v ) - cos(0, - dj)] 
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where ej = (I + l)/R are the energy levels for transverse vectors on S l x S 3 and 



V(6) = iJYQ - £ a, [log (l +e^-^ +w ) +log (l+e^ 



+Ar£v(e I .) 

(4.1) 
(4.2) 



where a/ = 2^-Z(Z + 1) and we have added the Lagrange multiplier JV to enforce the detP = 1 
constraint, i.e. 0, = 0. 

The equation of motion in terms of the Polyakov line eigenvalues z\ = e' 9i and the energies 
ejfl.e®, taking I lf =1 ^ / £ = &p( z ), becomes 



/=i 



z + e -p(p-ef) z + e -P(P+ef) 



z' + z £_2/(/ + 2)( z ' 2 -z 2 ) 



(4.3) 



z'-z £i{z'-zeP e r)(z!-ze-P e r) 

where V is the contour on which the eigenvalues of the Polyakov line lie (following [9, 10], but 
taking the contour to lie off the unit circle). In the confining regions of the p — T plane the contour 
is closed and the most general form of the eigenvalue density p(z), determined from the pole 
structure of the right hand side of eq. fOj, is given by 



p(z = a_ip_i + — + —j- + £ 
z ~- 



t>lPl | y blP-1 | y CI Pi 



(4.4) 



where the factors of p±\ indicate coefficients of terms with or without poles in c €, respectively. 
Solving the EOM and applying the constraint J ^f-p(z) = 1 gives the final form 



P(z) 
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(4.5) 
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where b tPl = -a, + 2t k =\ *(* + 2)fe/_*_ipi, (/ < L), and zj +1 = 2SJ =L+1 + l)e _/,e ' ' • The 
contour c € is given by z(j) and is obtained by inversion of p(z)dz = ids, which takes the form 



=zexp £ 
k=\ 



(-1) 
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z ~k e -kPn ^ b[pie kPe{ 
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x exp £ 

4=1 
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(4.6) 



In the low temperature limit this reproduces the results in [1]. In the deconfined phase the eigen- 
value density p (z) develops a gap as /J. is increased towards an energy level ej^ . To solve the EOM 
in this regime we can no longer use Cauchy's theorem in a straightforward way since the contour 
is not closed. Instead, we must define a resolvent 



fi>(z) 



1% 



dz' . ,,z + z' 



(4.7) 



and take the eigenvalues of the Polyakov line along a square root branch cut [9, 10]. The EOM is 
then obtained from the Plemelj formulae in terms of the average of the resolvent over the cut 



1 



zV'(z) = --[(o(z + s) + co(z-e)] , z e <T , 
and the density is obtained from the discontinuity of a>(z) across the cut 



(4.8) 



zp(z) = -[co(z + e)-co(z-e)] , z£tf. 



(4.9) 



Thus, observables can be computed from the resolvent by peeling the contour off the cut and 
collecting the surrounding poles using 



(4.10) 



It is now possible to calculate observables in both the confined and deconfined regions. The 
Polyakov lines can be obtained from &>±i = §^ ^§- i p{z)z ±1 in the confined regions, and from an 
expansion of the resolvent G)(z) = Tl T2£~ =1 z IFl ^±i in the deconfined regions. The quark 
number (normalized by A^ 2 ) is given by the Lagrange multiplier ,JV in the large N limit (as shown 
in [1]), which can also be calculated in both the confined and deconfined regions. In the confined 
regions it is found by solving the equation of motion to obtain the density p (z) and applying the 
normalization condition / ||p (z) = 1 to get an equation for jY . In the deconfined regions the EOM 
is solved for the resolvent o(z), and the det/ 5 = 1 constraint / ^^ft)(z)logz = gives an equation 
for JV . The full details and results for the calculation of both the Polyakov lines and the quark 
number are presented in [1] and resemble what is observed for N = 3, with the notable difference 
that the large Af results reveal more sharply defined transitions. As for the order of the transitions 
in the large Af theory, the normalized quark number JV is continuous as a function of the chemical 
potential. Also djVjd\i is continuous. But the third derivative of the potential, or d 2 jV /djl 2 , is 
discontinuous, indicating that the transitions are third order, of the Gross-Witten-Wadia type [9, 10]. 
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Figure 2: The phase diagram in the (fiR, TR) plane for QCD at large N and Nf for zero quark mass (Left) 
and large quark mass (Right). 
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Figure 3: Lattice results from N = 2 QCD for the quark number, Polyakov line, and quark number suscep- 
tibility as a function of the chemical potential, simulated on a 64 x 3 3 lattice with j3 g = ^= 24. 

The phase diagram in the jj.R - TR plane can be obtained from the equations for the density 



p(z) (eq. 4.5) and the contour z(s) (eq. 4.6) of the Polyakov line eigenvalues in the confining 
regions. The equations r(ju) for the lines of transition are obtained by calculating the z- values that 
give p (z) = 0, such that a gap is formed in the eigenvalue distribution, and then plugging these into 
the equation for z(s) such that the gap in the distribution is constrained to lie on the contour *6 '. 
The results for small and large quark mass are presented in Figure ||[ indicating that the series of 
confinement-deconfinement transitions is a low temperature feature which is delayed as a function 
of the chemical potential for non-zero quark mass, until that mass is (approximately) reached. 

5. Lattice results 

Lattice results are obtained from simulation of N = 2 QCD at non-zero chemical potential 
and are presented in more detail in [7]. The main result is reproduced here in Figure ^[ These 
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simulations were performed on a 64 x 3 3 lattice with a gauge coupling of j8 g = ^ = 24. The open 
and closed data points refer to different initial configurations and indicate the existence of multiple 
stable, or metastable, states. As seen for the theory on S l x S 3 , the quark number exhibits a level 
structure and the Polyakov lines show a deconfmement spike at each level transition. It is interesting 
to note the differences though. From the lattice, the curved shape of the steps in the quark number 
and matching behaviour in the Polyakov lines could be due to formulating the lattice theory on 
the torus, or it could be a result of working at stronger coupling. Another interesting feature, 
which is also observable perturbatively on S 1 x S 3 [7], is that the quark number susceptibility 
X q ~ TdjVjd\i follows the behavior of the Polyakov line, such that it also serves as an indicator 
of confinement-deconfinement transitions for \i ^ 0. 

6. Relationship between the quark number and Polyakov line 

Finally it is interesting to point out the relationship between the quark number and the Polyakov 



line which can be obtained in the large mass limit. From eq. ( |3.2| ), the quark number density goes 
to 

V3 2n 

in the large m limit, which agrees precisely with what is obtained in the lattice formulation of TV = 2 
QCD at non-zero density in [ll] 1 . It would be nice to explore this further. 

References 

[1] S. Hands, T. J. HoUowood and J. C. Myers, JHEP 1007 (2010) 086 [arXiv: 1003.5813 [hep-th]]. 
[2] P. de Forcrand, PoS LAT2009 (2009) 010 [arXiv: 1005.0539 [hep-lat]]. 

[3] A. Kurkela, P. Romatschke and A. Vuorinen, Phys. Rev. D 81 (2010) 105021 [arXiv:0912.1856 
[hep-ph]]. 

[4] G. Aarts and F. A. James, arXiv: 1009.5838 [hep-lat]. 

[5] O. Aharony, J. Marsano, S. Minwalla, K. Papadodimas and M. Van Raamsdonk, Adv. Theor. Math. 
Phys. 8 (2004) 603 [arXiv:hep-th/0310285]. 

[6] T. J. Hollowood and J. C. Myers, JHEP 0911 (2009) 008 [arXiv:0907.3665 [hep-th]]. 

[7] S. Hands, T. J. Hollowood and J. C. Myers, arXiv: 1010.0790 [hep-lat]. 

[8] B. Sundborg, Nucl. Phys. B 573 (2000) 349 [arXiv:hep-th/9908001]. 

[9] D. J. Gross and E. Witten, Phys. Rev. D 21 (1980) 446. 

[10] S. Wadia, 'A Study Of U(N) Lattice Gauge Theory In Two-Dimensions," EFI-79/44-CHICAGO. 

[11] K. Langfeld and G. Shin, Nucl. Phys. B 572 (2000) 266 [arXiv:hep-lat/9907006]. 



We thank Kurt Langfeld for pointing out this interesting result from the lattice. 



7 



